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Herschel-Bulkley 流体 的 
修正 多 松弛 参数 格子 玻 未 兹 曼 万 法 


Xf. WY, wS 
(南京 工业 大 学 机 械 与 动力 工程 学 院 ， 南 京 211800) 
摘要 : Herschel-Bulkley 流体 作为 非 牛顿 流体 中 重要 的 的 一 类 ， 已 广泛 应 用 于 工程 领域 。 为 了 改善 LBM 对 Herschel-Bulkley 流 
体 流 动 模拟 时 的 稳定 性 和 精度 ， 在 多 松弛 时 间 参 数 LBM(MRT LBM) 的 基础 上 提出 了 一 种 改进 方法 。 以 Poiseuille 流 为 例 ， 对 改 
进 算法 进行 了 详细 介绍 ， 对 剪 切 变 稀 型 和 前 切 增 稠 型 流体 分 别 进行 对 比 ， 并 比较 在 不 同 指数 情况 下 的 相对 误差 随 格子 数 增 大 的 
变化 趋势 ， 说 明了 此 改进 方法 的 有 效 性 ， 最 后 通过 在 水 泥 3D 打印 喷头 中 的 应 用 ， 与 多 松弛 时 间 参 数 的 LBM(MRT LBM) 的 模 
拟 结 果 进 行 对 比 ， 进 一 步 验证 了 此 改进 方法 的 可 行 性 。 
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A Modified Multiple-relaxation-time 
Lattice Boltzmann Method for Herschel-Bulkley Fluids 


WU Wei-Wei, HUANG Xiao-Diao, YUAN Hong 
(School of Mechanical and Power Engineering, Nanjing Tech University, Nanjing 211800, China) 
Abstract: As an important kind of non-Newtonian fluids, Herschel-Bulkley fluids have been applied into engineering 
widely. When the lattice Boltzmann method (LBM) was used in Herschel-Bulkley flow, a modified method based on 
multiple-relaxation-time lattice Boltzmann method was presented to improve the stability and accuracy of the 
simulation. Then Poiseuille flow was taken as an example to discuss the detailed process of the modified method. 
Also the variation trends of the relative errors with the lattice nodes in case of different power-law indexes were 
analyzed. The results showed that the modified method was efficient. Finally, the method was applied into the flow 
in the extruder of cement-3D printing, and the numerical results were compared with that of multiple-relaxation-time 


lattice Boltzmann method (MRT LBM), which further proved the feasibility of the method. 
Keywords: Herschel-Bulkley fluids; modified MRT LBM; Poiseuille flow; cement-3D prnting; MRT LBM 


0 引言 对 比 ， 结 果 验 证 了 此 方法 的 有 效 性 站，Malaspinas O 
"DENM "T 利用 局 部 方法 分 别 对 窜 律 流体 和 Carreau 流 体 进 行 了 

LBM 为 流体 流动 的 模拟 提供 了 除 有 限 元 力 法 之 。 详细 研究 ， 结 果 表明 在 一 定 的 指数 范围 内 ， 该 方法 
外 的 又 一 种 行 之 有 效 的 方法 ， 物 理 过 程 清晰 ， 计 算 。 有 较 好 的 稳定 性 ml Boyd J 等 采用 局 部 方法 对 备 律 流 
过 程 简 单 ， 已 被 泛 用 于 许多 复杂 流体 流动 中 > 人体 在 管道 中 的 流动 进行 了 分 析 ， 研 究 表明 该 方法 的 
Hr 对 非 和 牛顿 流体 流动 的 仿真 也 有 了 较 深 入 的 研 计算 精度 和 效率 均 高 于 Gabbanelli S 所 提出 的 截断 模 
完 。 在 分 析 非 牛 频 流体 流动 时 ， 松弛 时 间 会 随 着 流 型 [10], 2001 Æ, D'Humieres D 等 针对 三 维 方 腔 流 动 
体 当 地 粘度 的 实时 变化 而 变化 。 众 所 周知 ， 松 弛 时 提出 了 多 松弛 时 间 参 数 的 LBM， 并 证 明了 具有 上 比 普 


HAMRE KEMINEN, MABARA 解决 非 牛顿 流 体 的 仿真 问题 ，CHAI Zhenhua 等 利用 
精度 较 差 的 情况 “。 为 了 解决 这 个 问题 ， 相 关 学 者 都 MRT LBM 对 非 牛顿 流体 中 的 宕 律 流体 及 究 汉 流体 


提出 了 一 系列 的 方法 。Gabbanelli S EXI RIEUR E 进行 了 分 析 ， 结 果 表 明 此 方法 具有 很 好 的 稳定 性 和 
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Jit e fe Lk n SEAT f ERA. DAT T PERA E 
ERZO EFE; CHEN Songgui 等 采用 
MRT LBM 对 宾 汉 流体 在 平行 平板 中 的 三 维 流动 进 
fT Y EU. 解决 了 传统 LBM 对 非 牛 顿 流 体 仿真 时 的 
不 稳定 性 问题 叶 ，LI Qiuxiang 等 对 寡 律 流体 在 二 维 
方 腔 中 的 流动 进行 了 分 析 ， 结 果 表 明 MRT LBM 可 
以 很 好 地 应 用 于 对 非 牛 顿 流体 的 仿真 ， 而 雷诺 数 大 
小 和 军 律 指数 对 模拟 结果 有 较 大 的 影响 (1。 
虽然 对 非 牛顿 流体 的 研究 已 经 很 深入 ， 但 是 主 
要 还 是 集中 在 窜 律 流体 、 宾 汉 流 体 等 常见 的 非 牛顿 
流体 种 类 ， 对 于 不 常见 的 如 Herschel-Bulkley 流体 等 
研究 的 就 很 少 。 本 文 以 Herschel-Bulkley 流体 作为 
究 对 象 ， 基 于 MRT LBM 方法 ， 提 出 了 一 种 改进 方 
法 ， 可 以 更 有 效 地 改善 仿真 过 程 的 稳定 性 。 以 
Poiseuille 流 为 例 ， 详 细 介 绍 了 此 方法 ， 并 与 解析 解 
进行 了 对 比 ， 分 析 了 相对 误差 ， 最 后 将 此 方法 应 用 
于 对 水 泥浆 体 在 3D 打印 喷头 流动 中 的 分 析 ,， 对 该 方 
法 的 可 行 性 做 进一步 说 明 。 
1 MRT LBM 对 Herschel-Bulkley 流 体 的 
流动 模拟 
1.1 Herschel-Bulkley 流体 的 本 构 方程 
Herschel-Bulkley 流体 的 本 构 模 型 是 以 分 段 函数 
形式 呈现 的 ， 由 于 存在 初始 屈服 应 力 ， 当 小 于 初始 
屈服 应 力 时 ， 剪 切 率 为 0， 有 具体 表达 形式 为 : 


To+K7” ,lr|>7 
T,— 
7=0 ,人 <m 


(1) 


23 T E a Bt R AE ATT ERE, KR 
用 其 修正 本 构 方程 如 下 : 


7,71, [1 - exp(-my)] * Ky" (2) 


式 中 参数 m 主要 控制 应 力 的 增长 ， 来 避免 本 构 
方程 固有 的 非 连续 特性 ， 当 m 大 小 趋向 于 0 时 ， 即 
成 为 蝴 律 流体 的 本 构 方程 ， 当 m 趋向 于 无 穷 大 时 ， 
即 成 为 理想 的 Herschel-Bulkley 流体 ， 当 y=0 时 , 与 
原本 构 方 程 吻合 ， 在 一 般 情况 下 ，m 取 大 一 点 的 值 ， 
应 当 注 意 m 值 不 能 取 任 意 足 够 大 的 值 ， 否 则 会 带 来 
收敛 问题 ， 进 一 步 可 以 推导 得 到 显 观 粘性 系数 方程 


[16-18] 
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1.2 MRT LBM 流体 分 析 
MRT 格子 Boltzmann 7j ££ 7j: 


fr ejót,.t - ót) - f (r.t) - 


-MSM| f(r,)- f” G0] e 


Tg far DR f" rit) REEE Ar, fa (r,t), 


BPA, - Mf(r.r) , ü" (rn) - Mf" (r,t), ERE M E 


义 如 下 : 


1 1 1 1 1 1 1 1 1] 
4 -1 -1 -1 12 2 2 2 
4 23 .221 1 1 1 
0 1 0 -1 0 1 -1 -1 1 
M=I0 -2 0 2 0 1 -1 -1 1 
0 0 1 0 1 1 1 -1 -l 
0 0 -2 0 21] 1 -1 -I 
0 1 -1 1 -10 0 0 0 
0 0 0 0 0 1-1 1 -I 


S 为 与 松弛 相关 的 主 对 角 线 矩阵 : 


S = diag(5,,5,,55,54,54,55,55,55,5,) (5) 


ERAH, AEMET TEE, D 
度 和 动量 相关 的 量 so s3 ss IE 0, s; sgHX l/r 
这 里 的 tz 为 BGK 模型 里 的 松弛 时 间 ， 而 剩余 的 其 它 
参数 取 稍 大 于 1 的 数 ， 这 里 si=s=1.4，s4=s6c=1.2。 
MRT 模型 的 演变 过 程 由 两 个 过 程 构成 : 磁 撞 和 迁移 ， 
此 处 的 碰撞 步 与 BGK 模型 有 所 区 别 ， 按 下 述 方程 ; 
行 : 


广 (0D= fE, -M SER, D- “(rN)] (6) 


其 中 frANK E 2 BEAR ER E 
与 BGK 模型 的 相同 : 


fr * ejót,t - óty-f (r,t) (7) 


与 BGK 模型 相 比 ，MRT 模型 主要 包含 了 矩 空 
间 和 速度 空间 之 间 的 转换 过 程 。 利 用 
Chapman-Enskog 展开 ， 可 以 得 到 运动 粘度 : 


1 1 
u-pv-pc;—- -)àt (8) 
S 2 


也 可 以 得 到 密度 和 速度 的 宏观 量 定义 ， 与 BGK 模型 
相似 : 


p- 2, fr D. pu 2 ef erst) (9) 


运动 粘度 与 剪 切 速率 有 关 ， 传 统 数值 方法 求解 ， 
需要 建立 合适 的 差分 模型 , 而 LBM 可 以 避免 这 个 过 


T, fE MRT 模型 中 ， 应 变 率 张 量 比 BGK 模型 复杂 ， 
可 以 推导 得 到 02P9; 
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1.3 修正 MRTLBM 

根据 上 面 的 分 析 可 以 知道 , MRT LBM 应 用 于 非 
牛顿 流体 时 ， 主 要 受 当 地 松弛 时 间 的 影响 ， 也 可 看 
成 是 受 粘 度 的 影响 ， 而 粘度 通过 应 力 张 量 与 应 变 率 
相关 联 。 对 于 Herschel-Bulkley 流体 ， 当 n<1 时 ， 对 
应 剪 切 变 稀 型 流体 ， 当 n=1 时 ， 则 为 宾 汉 流体 ， 当 
n>1 时 ， 对 应 剪 切 增 稠 型 流体 。 在 对 含有 指数 项 的 非 
牛顿 流体 直接 应 用 MRT LBM 分 析 时 会 出 现 稳定 性 
和 精度 问题 ， 当 剪 切 率 趋 近 于 0 时 ， 对 于 n<1 的 剪 
切 变 稀 型 流体 ， 粘 度 会 趋向 于 无 穷 大 ， 导 致 发 散 ; 
对 于 n>1 的 剪 切 增 稠 型 流体 ， 粘 度 会 趋向 于 0. 在 过 
去 的 学 者 研究 过 程 中 ， 并 未 针对 Herschel-Bulkley jfi 
体 的 分 析 提 及 相关 的 解决 方法 。 
事实 上 ， 大 部 分 的 非 牛顿 流体 仅仅 在 某 一 个 剪 
切 率 范围 内 呈现 出 非 牛顿 流体 的 特性 ， 在 此 范围 外 
时 ， 通 常 表现 出 类 似 牛 顿 流体 的 特性 ， 即 粘度 是 一 
个 定常 数目 ， 所 以 基于 这 种 现象 ， 针 对 Herschel- 
Bulkley 流体 ， 提 出 了 一 种 修正 MRT LBM 方法 ， 这 
种 方法 以 分 段 函 数 的 形式 呈现 。 由 于 
Herschel-Bulkley 流体 存在 屈服 应 力 ， 为 了 保证 剪 切 
率 在 有 效 范 围 外 时 ， 粘 度 是 一 个 常数 ， 这 里 用 宾 汉 
流体 的 本 构 模型 来 进行 代替 。 


t, / Yo bud" PASA 
W= To AKI? hiia (11) 
Tol Jo FUS” Y. 


MRT LBM 可 以 在 动力 粘度 处 于 某 一 范围 时 , 进 
行 精确 、 稳 定 地 模拟 分 析 ， 而 当 松 弛 时 间接 近 于 1/2 
时 ， 会 出 现 不 稳定 的 情况 ， 引 起 发 散 ， 此 时 对 应 的 
粘度 接近 于 0, v<0.001; 对 应 地 ， 当 松弛 时 间 大 于 1 
时 ，MRT LBM 的 精度 会 降低 ， 此 时 对 应 的 粘度 相对 
较 大 ,之 1/6。 因 此 对 式 (11) 设 置 粘度 的 上 极限 值 和 下 
极限 值 wwan=0.001 和 vwax=0.16。 很 明显 , 对 于 剪 切 变 
稀 型 流体 ， 粘 度 最 大 值 对 应 剪 切 率 接 近 于 0 的 时 候 ， 
相应 地 ， 对 于 剪 切 增 笛 型 流体 ， 粘 度 最 小 值 对 应 剪 
切 率 接近 于 0 的 时 候 。 
2 平行 平板 之 间 的 流动 
由 于 泊 肃 叶 流 存在 理论 解 ， 常 被 用 来 验证 稳定 
性 及 精度 问题 。 这 里 以 泊 肃 叶 流 为 例 ， 介 绍 这 种 修 
正 模型 的 具体 过 程 。 两 平行 平板 沿 y 方向 的 距离 为 
五, 在 x 方向 存在 一 个 恒定 值 的 压力 梯度 VP， 可 以 
得 到 沿 y 方向 的 速度 分 量 wu 的 分 布 情况 。 根 据 对 称 
性 ， 将 上方 向 的 距离 分 成 [-2V2, 0] 和 [0, 8/2] 两 部 分 ， 
根据 对 称 性 ， 只 需 对 [0, /2] 范 围 内 分 析 ， 就 可 以 得 
到 久 沿 整个 y 方 向 的 分 布 情况 。 在 [0, HW/2] 区 间 内 ， 
可 以 分 出 四 个 不 同 的 区 域 A、B、C 和 D， 其 中 区 域 
A 对 应 的 是 靠近 壁面 的 高 剪 切 率 范 围 , 此 时 剪 切 率 超 
过 了 yp,， 以 yi 作为 分 界 值 ， 在 此 区 域内 ， 将 流体 看 
成 宾 汉 流体 ; 区 域 B 对 应 的 是 Herschel-Bulkley 流体 
的 有 效 区 域 ， 以 y, 作 为 临界 分 界 值 ， 区 域 C 对 应 的 
是 低 剪 切 率 范围 ， 此 时 的 流体 看 成 是 宾 汉 流体 ， 区 
B D 对 应 的 是 应 力 小 于 屈服 应 力 的 阶段 ， 此 时 剪 切 
率 为 0, 以 y 作 为 分 界 值 .结合 泊 肃 叶 流 的 速度 推导 ， 
可 以 得 到 该 修正 方法 的 速度 分 布 的 理论 函数 。 
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上 式 中 的 常数 项 wl、o 和 o5 RARU F: 
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2d (n To El NE? ) þa 
avp | |2 oPiæ) V dp/ : 
Herschel-Bulkley 流体 相对 应 的 延续 ， 垂 直 的 虚线 代 
表 了 不 同 区 域 的 分 界 点 。 
-5 
x10 


将 上 述 速度 分 布 的 函数 求 导 后 ， 根 据 应 力 相 等 
条 件 可 以 得 到 对 应 的 分 界 值 分 别 为 : 


—. fo 
OP / Ox 


, y l/(n-1) 
PVo 
^ zs K | d 04) 


n NVO) 
" 1 (pva) : 
^ — 8P/éx K d 


其 中 6P/6x=VP ， 根 据 表达 式 中 物理 量 的 数值 
大 小 的 不 同 ， 可 以 得 到 不 同 的 区 域 划分 。 
根据 式 (14)， 这 里 主要 考虑 四 个 因素 对 最 后 仿真 
结果 的 影响 :初始 屈服 应 力 ww， 压 力 梯度 VP， 窜 律 
指数 n 以 及 格子 数 N。 观察 式 (14) 后 发 现 ， 为 了 保证 


NM-D) 
yy 或 (和 )ys 的 存在 ， matr «(Eo E 


i 703 y l/(n-1) 
(和 ) To < sfer) o 


为 了 考察 使 用 此 截断 模型 时 ， 不 同 的 指数 n 的 
有 效 性 ， 这 里 取 的 范围 为 0.2~5， 这 里 窜 律 指数 n 
取 0.5 和 1.5， 分 别 对 应 剪 切 变 稀 型 流体 和 前 切 增 秽 
型 流体 ， 这 样 可 以 更 清楚 地 表现 初始 屈服 应 力 对 仿 
真 结果 的 影响 ,。 当 n 取 0.5 时 ,压力 梯度 VP=-1.5x10， 
流体 的 粘度 系数 K=0.001， 平 板 的 高 度 10， 在 模 
拟 时 ， 格 子 数 取 为 200， 流 体 的 密度 p=1， 运 动 粘度 
的 极限 y=0.001, v,,4,70.16, 这 里 n<1, MÆ v0.16, 
yw=0.001， 初 始 届 服 应 力 wo=3.0x10*， 其 数值 结果 和 
仿真 结果 对 比 图 如 图 1 所 示 ， 图 中 实 线 是 根据 式 (12) 
的 解析 解 描绘 的 速度 分 布 情况 ,圆圈 是 通过 LBM 仿 
真 获 得 的 ， 虚 线 部 分 分 别 是 宾 汉 流体 和 


yr = 


-5 -yl-yt 0 yvy 5 


y 
图 1 n-0.5 时 的 速度 分 布 
Fig. 1 Velocity distribution for n=0.5 


当 究 律 指数 n 取 1.5 时 ， 压 力 梯度 VP=-2x10”， 
流体 的 粘度 系数 K=0.01, 平板 的 高 度 H=10, 在 模拟 
时 ， 格 子 数 取 为 200， 流 体 的 密度 p=1， 运 动 粘度 的 
极限 y=0.001,; v,,4,70.16, 这 里 n>1, MÆ ve=0.001， 
w=0.16， 初 始 届 服 应 力 w=4.0x10*， 其 数值 结果 和 
仿真 结果 对 比 图 如 图 2 所 示 ， 图 中 实 线 是 根据 式 (12) 
的 解析 解 描绘 的 速度 分 布 情况 ,圆圈 是 通过 LBM 仿 
真 获得 的 ， 虚 线 部 分 分 别 是 宾 汉 流体 和 
Herschel-Bulkley 流体 相对 应 的 延续 ， 垂 直 的 虚线 代 
表 了 不 同 区 域 的 分 界 点 。 

观察 图 1 和 图 2 可 以 发 现 ， 无 论 初始 屈服 应 力 
如 何 变化 ， 数 值 解 与 解析 解 具 有 较 好 的 重合 度 。 初 
始 届 服 应 力 w 与 压力 梯度 YP 决定 了 届 服 阶段 即 区 域 
D 的 范围 大 小 。 
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图 2 n-1.5 时 的 速度 分 布 对 比 
Fig. 2 Velocity distribution for n=1.5 


这 里 对 于 四 种 指数 不 同 的 流体 ， 均 进行 了 一 系 
列 的 仿真 ， 比 较 MRT LBM 仿真 与 解析 解 之 间 的 相 
对 误差 ,其 中 相对 误差 通过 > (4-wv /v A 
计算 获得 。 由 于 LBM 仿真 结果 的 精度 与 粘度 有 关 ， 
为 了 保证 只 讨论 相对 误差 与 格子 数 之 间 的 关系 ， 在 
进行 单位 转换 的 时 ， 模 型 中 的 粘度 v'—w(0x 0n, 4 
6x t, HK LBM 的 格子 数 N 时 ， 只 需要 同时 按 
比例 缩放 5x FU ôt, 即 可 保证 粘度 不 变 ,由 于 5xrxLN ， 
根据 前 面 的 关系 可 知 ，5tw1/N?， 这 里 入 逐渐 从 10 
增 大 到 200, 在 指数 和 压力 梯度 不 同 的 情况 下 ， 相 对 
误差 如 图 3 所 示 。 


1 
1 2 
10 10 


格子 数 
3 相对 误差 对 比 图 


Fig. 3 Relative Error analysis 

图 中 实 线 代表 了 和 斜率 为 -1 的 直线 ， 图 例 中 a 表 
示 低 剪 切 率 区 域 范围 C 较 小 的 情况 ,b 表示 低 剪 切 率 
区 域 C. 范围 较 大 的 情况 , 对 比 发 现 , a RI b 两 种 情况 
下 产生 的 相对 误差 基本 相同 ， 且 随 着 格子 数 的 增加 ， 
相对 误差 不 断 地 下 降 ， 如 果 进 一 步 增加 格子 数量 ， 
可 以 控制 在 0.01% 以 内 ， 从 图 中 可 以 进一步 看 出 ， 当 
指数 取 值 不 同时 ， 相 对 误差 的 差异 性 较 明显 。 
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3 在 水 泥 3D 打印 喷头 内 的 仿真 应 用 


在 进行 水 泥 3D 打印 时 ， 需 要 利用 螺杆 挤 出 方式 
来 实现 水 泥浆 体 的 挤 出 成 型 ， 因 此 了 解 水 泥浆 体 在 
喷头 内 的 流动 显得 尤为 重要 。 水 泥浆 体 在 一 定 的 7 
灰 比 下 呈现 出 明显 的 Herschel-Bulkley 流体 特性 ， 当 
室温 为 20 'C， 水 灰 比 为 0.5， 水 化 时 间 为 5 min 时 ， 
浆 体 呈现 出 典型 的 Herschel-Bulkley 流 变 模式 ， 流 变 
学 方程 描述 为 0 


[3.899+1.1037°® ,lr,|> 3.899 
7=0 J[r,| «3.899 


(15) 


将 上 式 转 换 成 式 (11) 的 形式 , 式 中 参数 m 主要 控 
制 应 力 的 增长 ， 这 里 取 m=500， 则 有 


7—3.899[1— exp(—m7)]+1.1037” (16) 


水 泥浆 体 在 图 4 的 螺 模 区 域 流动 ， 其 中 螺杆 以 
恒定 的 转速 N 进行 旋转 ， 假 设 螺 杆 静 止 ， 则 可 认为 
是 机 简 以 转速 N 进行 反 向 旋转 ， 将 整个 结构 进行 展 
开 ， 即 可 形成 一 个 截面 为 矩形 的 腔 情 ， 此 时 忽略 螺 
棱 的 轻微 倾斜 请 ， 展 开 示意 图 如 图 5 所 示 ， 对 图 4 
水 泥浆 体 在 螺杆 中 的 流动 分 析 即 转换 成 水 泥浆 体 在 
图 5 中 长 方 体型 腔 体 中 的 流动 分 析 。 
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图 4 喷头 结构 图 
Fig. 4 Extruder Structure 
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图 5 喷头 展开 示意 图 
Fig. 5 Expansion Structure 


根据 喷头 的 展开 示意 图 ， 这 里 取 中 心 位 置 的 截 
IBi(Y-Z)ceue 采用 D2Q9 模型 ， 其 中 只 对 上 表面 设置 
速度 ， 速 度 方向 治 y-h 处 的 Z 方向 ， 其 他 三 个 面 速 
度 均 为 0， 忽略 壁面 滑 移 的 影响 ,截面 的 实际 尺寸 为 
W-16 mm, /-6 mm， 结 合 喷头 扫描 速度 的 需求 以 及 
水 泥浆 体 的 凝固 时 间 ， 浆 体 的 挤 出 速度 相对 较 低 ， 
螺杆 的 旋转 速度 取 30 r/min 进行 ， 螺 杆 的 导 程 角 
0=20 ， 利 用 相似 准则 ， 取 雷诺 数 Re 为 关键 准则 数 
进行 量 纲 转换 ， 可 以 得 到 仿真 模型 中 的 各 参数 ， 再 
j 式 (13) 对 粘度 进行 修正 ， 即 可 保证 稳定 性 和 精 


Jt 


采 


ER 


度 ， 得 出 流 线 图 如 图 6 所 示 。 流 体 的 流动 以 环流 的 运用 修正 MRT LBM, 采取 不 同 的 格子 数 进行 模 
形式 呈现 ， 环 流 的 中 心 在 (0.5 刺 ，0.7 入 附近 ， 在 左下 拟 ， 并 与 MRT LBM 进行 比较 ， 得 到 相应 的 速度 场 
角 和 右 下 角 不 存在 明显 的 流体 流动 。 分 布 ， 如 图 7 所 示 ， 其 中 修正 MRT LBM 格子 数 分 
别 取 288x108 和 160x60，MRT LBM 格子 数 取 
288x108 与 修正 MRT LBM 进行 比较 ,图 7 (a) 是 取 螺 
E si HEEE 8 mm 时 , WENE u 与 螺 覃 深度 之 间 的 关 
2 34 AE, 图 7 ORIRE 5.4 mm 时 ， 速 度 分 量 
FINI u 与 螺 槽 宽度 之 间 的 关系 图 , 图 7 (9) 是 取 螺 槽 宽度 为 
14.4 mm 时 ， 速 度 分量 v 与 螺 模 深度 之 间 的 关系 图 ， 
| 图 7(d) 是 取 螺 槽 深度 为 3 mm 时 , 速度 分 量 v 与 螺 模 
0 d E D ik 宽度 之 间 的 关系 图 。 可 以 看 出 该 修正 方法 可 以 有 效 
EB 对 速度 场 的 分 布 进行 模拟 ， 采 取 不 同 的 格子 数 时 ， 
Pe Seanin 虽然 结果 存在 一 定 误差 ， 但 是 基本 一 致 。 
6 5 
3 4 
E^ gs 
&» ——  MRTLBM 288*108 | 8 1 — MRTLBNOSS* 108 
— 288*108 ----- 28S *108 
loW gus 160*60 OV se 160*60 
0 : ; ; -1 : : 
-2 0 2 4 6 8 0 4 8 12 16 
u/mm-:s-1 Width/mm 
图 7(a) 速度 wu- 深度 有 h 图 7(b) 速度 u- ERE W 
Fig. 7(a) Velocity u-Depth A Fig. 7(b) Velocity u-Width W 
: 9 
JD 1 
BA NR. E 
3 al ume E 0 
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-1.5 -1 -0.5 0 0 4 8 12 16 
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图 7(c) 速度 vw 深度 有 h 图 7(d) 速度 宽度 W 
Fig. 7(c) Velocity v-Depth A Fig. 7(d) Velocity v-Width W 
图 7 速度 分 布 图 


Fig. 7 Velocity distributions 


4 结论 


针对 广义 牛顿 流体 中 的 Herschel-Bulkley 流体 ， 
提出 了 一 种 修正 MRT LBM 方法 ， 可 以 有 效 改 善 在 


模拟 时 易 出 现 发 散 的 情况 。 
(1) 通过 Poiseuille 流 对 此 方法 进行 了 详细 


论 


xs, H 
子 数 的 增加 ， 
且 剪 切 变 稀 型 和 剪 切 增 稠 型 流体 呈现 相同 的 变化 趋 
势 ， 还 进一步 证 明了 低 剪 切 率 
对 误差 几乎 没有 影响 ， 而 指数 n 


N 


导 了 理论 解 公式 ， 计 算 了 相对 误差 ， 随 着 格 


相对 误差 不 断 减 小 , 能 达到 0.05% 左 右 ， 


区 域 范 围 的 大 小 对 相 
的 取 值 不 同 对 其 有 


较 大 影响 。 
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C20 将 此 方法 应 用 于 对 水 泥浆 体 在 3D 打印 喷 


头 中 流动 情况 的 分 析 ， 并 采用 MRT LBM 方法 与 文 


中 的 方法 进行 了 对 比 ， 发 现 其 速度 场 分 布 基本 一 致 ， 


步 证 明了 此 改进 方法 的 有 效 性 。 
(3) 对 水 泥 3D 打印 喷头 中 的 流动 分 析 可 知 ， 


受 浆 体 粘性 和 喷头 结构 的 影响 ， 浆 体 在 喷头 横 截 面 
的 流动 呈现 环流 特征 ， 环 流 中 心 大 致 在 (0.5WW，0.7 有 h) 


处 ， 


可 以 根据 流 线 图 来 判断 浆 体 的 流动 区 域 。 从 流 


动 分 析 可 以 看 出 ， 在 腔 体 的 下 方 左右 两 角 处 不 存在 
明显 的 流体 流动 ， 而 为 了 保证 流体 的 流动 速度 ， 一 
方面 可 以 适当 增 大 螺杆 的 旋转 速度 ， 另 一 方面 可 以 
适当 增 大 螺 模 的 宽度 。 
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